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The Flory theory for a single polymer chain is derived as the lowest order of a cumulant 
expansion. In this approach, the full original Flory free energy ( including the logarithmic term ), 
^ I is recovered. The prefactors of the elastic and repulsive energy are calculated from the microscopic 

. parameters. The method can be applied to other types of monomer-monomer interactions, and the 

■ case of a single chain in a bad solvent is discussed . The method is easily generalized to many chain 

' systems (polymers in solutions), yielding the usual crossovers with chain concentration. Finally, 

, this method is suitable for a systematic expansion around the Flory theory. The corrections to 



Flory theory consist of extensive terms ( proportional to the number N of monomers ) and powers 
of A^2-!^d xhese last terms diverge in the thermodynamic limit, but less rapidly than the usual 



' Fixman expansion in '^1'^ . 
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LA THEORIE DE FLORY REVISITEE 



RESUME 



La theorie de Flory pour une chaine polymerique est obtenue comme I'ordre dominant d'un 
developpement en cumulants. Dans cette approche, I'energie libre originale de Flory ( y compris 
le terme logarithmique ) est obtenue. Les prefacteurs des energies libres elastique ct repulsive sont 
derives a partir des paramctres microscopiques. La methode pent etre appliquee a d'autres types 
d'interactions entre monomeres, et on discute le cas d'une chaine en mauvais solvant. La methode 
peut etre gencralisee au cas de plusieurs chaincs ( solutions de polymeres), et on en deduit les 
changcmcnts de comportcmcnt en fonction de la concentration en chaines. Finalement, la methode 
permet un developpement systematique autour de la theorie de Flory. Les corrections a la theorie 
de Flory comportent des termes extensifs ( proportionnels au nombre N de monomeres ) ct des 
puissances dc jV^"'^'^. Ces tcrmcs divergent a la limite thermodynamique, mais moins vite que le 
developpement de Fixman, en puissances de A^2-ci/2 
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I. INTRODUCTION 

The size of a polymer in a solvent is characterized by an exponent u, which relates the radius 
of gyration R of the chain to the degree of polymerization N (number of monomers) through the 
formula^^'^'^) 

R^N'' (1) 

In an ideal O solvent, chains are Brownian with i/ = 1/2. In a good solvent, chains are swollen, 
and the index u is larger than 1/2. 

Flory*^^-* has devised a simple argument to compute the exponent u of a chain in a good 
solvent, which gives amazingly good agreement for all dimensions d. The argument, as presented in 
ref.(2) goes as follows: consider a chain of monomers of length a, with local repulsive interaction 
characterized by an excluded volume parameter v. The swollen chain will acquire a radius of gyration 
R; the monomer concentration is c = N/R'^, and the repulsive energy is thus 



Er.^ = R'^ 

1 iV2 



^rep 



(2a) 



2^^ 



The elastic energy of the chain is given by: 



up to a multiplicative constant ( T is the temperature). 
The total free energy is: 

i?2 t,iV2 



and minimization with respect to R yields the celebrated Flory exponent: 

which is exact for d = 1,2, A, and almost exact for d = 3 (see ref.(3); up = 0.6 whereas best 
numerical estimates give u = 0.588). 

Note that this simplified derivation of the Flory theory misses a logarithmic term in the free 
energy, which is present in the original derivation of Flory (ref.(l) ). In any case, the derivations 
of the Flory exponent are quite empirical, and any attempt to go beyond it, by improving on any 
of the two terms of (2a) or (2b) turns out to ruin the argument (see ref.(2) for a discussion of that 
point). Moreover, the free energy increases like N instead of increasing like iV as it should in 
the asymptotic limit of large N (the free energy is not extensive) . 

Des Cloizeaux has proposed a Gaussian variational procedure to tackle this problem '^^^ but 
the outcome was i^^ = |, by far too large a result. Also, Edwards has proposed a self-consistent 
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method^^) which yielded the Flory formula, but, as noted by des Cloizeaux^^) , this approximation 
yields unrealistic chain configurations. 

In this paper, we propose a cumulant expansion, which to lowest order generates a Flory type 
free energy from the microscopic model, and which can be easily generalized to the study of more 
complex systems (polyelectrolytes, chains in bad solvents, membranes and interfaces, solutions of 
polymers, etc.). This method allows for a systematic expansion around Flory theory. 

II. ONE CHAIN SYSTEM 

We start from the Edwards continuous representation of a chain: 

Z = [ Vr{s) dl^f ds f(s) ) exp ( - J / ds P - ^ [ dsds'v (f(s) - f{s'))] (5) 

Jf'(N)=f{0) \^ Jo J \ ^ Jo ^ Jo J 

where V [r — r ') is the interaction between a monomer at point r and a monomer at r ' and units 
have been chosen so that the persistence length a is equal to 1. 

The first (5-funetion is used to constrain the center of mass of the chain at the origin. Further- 
more, to simplify calculations, we have assumed that the chain is closed, i.e. ri^N) = r(0). As will 
be seen later, this constraint is unessential, as far as the method is concerned. 

Performing a Gaussian transform on (5), and denoting by v~^{r) the inverse of the kernel 
u (r) , we have : 

Z = J ^0(^)exp^— -y dfdf '(j) (f) {f — f ') (j) (f (6a) 



where 

Zrh 



[ Vf{s) dl^f ds r(s) ) exp f - J / ds P - i [ ds (/)(r(s)) J (66) 

7r(Ar)=f(0) Jo / \ "'O "'O / 



Using the identity: 
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y dQ <5 - 1 ds f2(s) j (7a) 
where Q is the square of the radius of gyration, we can rewrite (6b) as: 

Z^= dQ Z^{Q) (76) 
Jo 



where 

MQ) = [ 1^r{s) 5 (1 r ds f(s) ) sIq-^ r ds r\s) 

Jf'{N)=r{0) Y Jo J \ -'^ Jo 

xexp(— -y ds P — i J ds^(r(s))j 



(7c) 



Defining (...) as the average with respect to the partition function Zo{Q) = Z^=o{Q), we can 
write , , 



^ d5 0(r(s))j 
= ^exp ^— i J df (j) (r) p (r)^ 



(8) 



where p (r) = ds 6 {f— f(s)) . 

We use the cumulant expansion: 



<e-^> = exp (-{A) + i «A^> - {Af) - i {{A^) - 3 {A^) {A) + 2{Af) + 



(9) 



To lowest order, we have: 



Z^iQ) ~ ZoiQ) exp (^-i j df4> (f) pQ (f) j (10a) 



where 
Z, 



(lOi.) 

and 



= 7^4^/ ds /" 1^r(s)6l^f ds r(s)] 5 I Q - ^ [ ds r^is)] 



xexp ( ~2 y ds ^(s) ) (r — r(s)) 



Zo{Q) and pQ{r) are respectively the partition function and monomer concentration for a Brownian 
chain with center of mass constrained at the origin and radius of gyration square constrained to Q. 
Finally, performing the {0} Gaussian integral of (6a), we obtain: 

Zc^Zf = J^ dQ ZoiQ) exp (^-^J drdr'pQir)vir-r')pQir')^ (H) 

The F index in (11) stands for Flory, since, as we shall see, Zp is just the Flory partition function. 

A straightforward calculation of (10c) and (lOd) yields ( see the appendix ), for Q/N — > oo 
(swollen chain) : 



and 

J2 



In the thermodynamic hmit, N — >■ +00, and (11) can be evaluated by applying the saddle- 
point method on Q. Assuming a standard contact interaction: 



V (r) = V 6 {r) (13) 
where v is the excluded volume parameter, the function to minimize is: 



Nj N 2 \ATr J Q'i/s 

where f3 is the inverse temperature, and the radius of gyration Rq = \fQ satisfies the minimization 
equation: 

Rg N 2 \Atx) ^ ' 

We recognize in equations (14) the Flory free energy, with calculated prefactors, including the 
logarithmic term, present in the original Flory theory. The role of this term is just to ensure a 
crossover from a Brownian regime in d > 4 to a swollen regime in d < 4. Indeed, to solve (14b), we 
must balance two terms of opposite sign, and check that the third one is negligeable with respect 
to them. 

Balancing the 2"^*^ and ?>^'^ term in (14b) yields the Flory result, = The consistency 
requirement: <^ is satisfied provided that d < 4. 

Balancing the 1*^' and 2"*^ term in (14b) yields the Brownian exponent u = 1/2, and the 
consistency check ^^-^ <C is satisfied for d > 4. 

Thus the logarithmic term indeed, enforces the crossover between the 2 regimes. 

The method can be applied to any type of two-body interaction v (r) and denoting by v{k) its 
Fourier transform, the Flory free energy becomes: 




(iFp = -{d - 1) In ^^j + 27r^^ - I y dkv{k) exp ( -Q- ] (14c) 

Application of this method to the case of polyelectrolytes, where the two-body interaction is the 
Coulomb potential v (r) = ^(d_2) > yields the well-known Flory exponent ^'^'^ u = 2 ■ This is known 
to be incorrect for dimensions larger than 3 , where the exact exponent is v = ( see rcf(7) ). 

Finally, let us show how this method works for a chain in a poor solvent. In that case, in 
addition to the second virial coefficient v which is attractive, one has to introduce the third virial 
coefficient w which is repulsive. The partition function of the chain reads: 



/ ^r{s)si^ rdsf{s) 

Ji'(N)=r(0) \ Jo 



li'{N)=f{0) 
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dsP + ^ dsds'S{r{s)-r{s'))-^ / dsds'ds"6 {r{s) - r{s')) S {f{s) - r{s")) 
^ Jo ^ Jo ^ Jo ^ 

(15a)' 

In the present case, the alternative to the Gaussian transform is to constrain the monomer density 
p(r) as a new integration variable through the identity: 

l = y I>(^(r)I>p(r) exp J dr(l){r)p{r) — i J ds (r(s))^ (156) 

Inserting identity (15b) in equation (15a) yields: 

Z = J V^{r)Vp{r)e^p (^i J dr^{f)p{r) + ^ j drp^{r) J drp^{f)^ (15c) 

where is defined by eq.(6b). Using identity (7a) in the above equation together with the lowest 
order cumulant ( see eq.(lOa)) and integrating over the fields (p and p yields the expression: 

Z ^ y dQZoiQ) exp (+^ J dfplir) -^J dfpl^{r)^ (16) 

where Zo{Q) and pQ (r) are defined in eq.(10 b,c). 

In the case v < 0, we are back to the previous situation, and the chain is swollen. The free 
energy reads: 

/3^F = -(c/-l)ln^ +2vr2^-- — + (17a) 



J N 2\AnJ 6 V187ry 

At the point, v = 0, and the free energy is still given by the above expression with v = . 
This regime corresponds to an exponent u = In d = 3, as is well known, this yields a Brownian 
exponent u = 1/2. Note that , in that case, since the free energy of equation (17a) is finite, one 
cannot evaluate the integral (16) by the saddle-point method. In that case, quadratic corrections 
must be included to compute the free energy. 

In the collapsed regime, Q/N — > and Zo{Q) is no more given by (12a), although the 
asymptotic form of pq (r) is given by (12b). We show in the appendix that the correct form for 
Zq{Q) in dimension d = 3 is: 

3 f N\^^'^ N 



yielding a different Flory free energy: 

^ ^~ ''U6^27r)Uy J^8Q 2^4^; Q^/^ + 6 
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In this regime, we recover the usual exponent = |- 

The method described above, can be applied to any type of interaction v (r) , and can be easily 
generalized to the case of membranes or interfaces 

III. PAIR CORRELATION FUNCTION 

At this level of approximation, it is easy to compute the corresponding pair correlation function 
g{f) . Using the following definition: 



= j^J^ ds ds' {5{f{s))5{f-f{s'))) (18a) 



or 



m = J^J^ ds' (exp(zg(f(s) - f(s')))) (186) 

in Fourier space. Using the lowest order cumulant expansion (10a), we show in the appendix that 
the pair correlation function is given by: 



d-l 



xexp(-:^(l + 2^-(^-«)sm9)) (19o) 

where the square radius of gyration Q is given by its saddle-point value as computed in the previous 
section, and J^-i denotes the Bessel function of integer order (c? — 1) . 
At small q, i.e. large distances (fQ « 1 , the function reduces to: 

g{q)=N — exp(-Q(l - cos e)f/d) (1%) 
whereas at smaller distances g^Q > 1 , a simple scaling argument shows that: 

in three dimensions. This result agrees with the prediction of Edwards^^-*. 

IV. MANY CHAINS SYSTEMS 

Consider now a system of M polymer chains of length N (polydispersity effects can be trivially 
included). The partition function of the system reads: 

» N / AT ^.N pN \ 

Z= \{Vr,{s)e^v\-^Y. - - ^ / dsds'5 {n{s) - f,{s))\ (20) 



Introducing a center of mass coordinate Ri and a square radius of gyration Qi for each chain, 
the same procedure as in II) yields the fohowing expression: 



=J llmdQi exp {(d - 1) In - 2n'^ 



d 



27r {Qi + Qj 



d/2 



cxp 



V 



2 Qi + Qj 



(21a) 



Each polymer chain is represented by a coordinate Ri (center of mass) and a radius of gyration 
square Qj,. The determination of tlie size Qi will be done by a saddle-point method, and it is natural 
(for monodisperse systems) to assume that all sizes arc equal: Qi = Q. We obtain: 

2Q' 



Zp = exp Af 



{d - 1) In 



Q 

N 



27r- 



N 

z=l \ ^ " / 1,3 = 1 



Rj — R 




(216) 



Thus, the partition function contains a Flory-Iike term, characteristic of each individual chain 
properties, and a partition function representing a liquid of interacting chains with a smoothed 
interaction: 



V [r^ = vN^T 



d 



47rQ 



d/2 



exp 



d ^2 



(22) 



The range of the interaction is the size of the chains \/2Q. Once Zp is calculated, the radius 
of gyration Rq is to be determined by minimizing the total free energy of the system with respect 
to Q. 

Defining the monomer concentration c = NN ( where Q is the volume ) , and the overlap 
threshold concentration c* = ~ jSf^'^-^d.) ( where Rg = \fQ is the radius of gyration of a chain, 
and V is the Flory exponent v = ) , the various regimes of chain concentration can easily be 

recovered. 

i) Dilute regime c « c*: it is the gas type regime, where the density of chains is such that 
the typical interchain distance is much larger than the chain size Rq . In that case, we have a 
gas of weakly interacting single chains. The partition function can be evaluated by using the virial 
expansion , and the osmotic pressure reads: 



(23a) 



where c/N = Af/Q is the polymer concentration, and A2 is the second virial coefficient calculated 
with (22): 

A2 = -^JdR (exp{-pV{R)) - 1) (236) 



A simple evaluation of this integral yields: 



A2C^R%{\nNY/^ (23c) 

Note the logarithmic term in the virial coefficient ^2, absent from standard calculations. The total 
free energy per polymer reads: 

F/M = Ff + C/NA2T {23d) 

where Fp is the single chain Flory free energy as given by (14a). Minimization of (23d) with respect 
to Q yields the usual Flory exponent. 

iii) Melt c*<<c2±l/a'^( where a is the monomer size ) : when interchain distances become 
smaller than the radius Rg, the chains overlap strongly, and the monomer concentration is constant, 
with small fluctuations. Using a Gaussian transform, equation (21b) can be written as: 



Zp = exp ( Af 



X J V(t){R) exp(^-^ J^^J d-^>(-R)^"^(^-^O'^(-^O+-^My"d-Rexp(i/30(i?)))^ 

(24a) 

where V-'^{R-R') is the inverse kernel oi V{R — R'). 

This functional integral can be evaluated by the saddle-point method. Assuming a uniform 
monomer concentration is equivalent to assume a constant field (f>{R) = (po- A. simple calculation 
shows that the saddle-point is given by: 

^(/)o = iNcv'^ (246) 

so that the free energy becomes: 

/Jft . -A^ ((<i - 1) ,„ I - 2.^1 + 1„ fi) + ^-ml ,24c) 

Minimization of (19c) with respect to Q yields the standard Brownian exponent u = 1/2 , and the 
osmotic pressure is given by: 11/? = f 

The interchain interactions screen out completely, and the chains become Brownian, as is well 
known(^°). A systematic expansion around the mean-field ^0 can easily be done, and the lowest 
order turns out to be identical to the usual RPA . 

ii) Semi-dilute regime c* « c « : it is the liquid type regime, which requires more 

involved calculations of the liquid-like partition function '^^^ . A careful analysis of the generic term 
of the virial expansion of (16b) leads to the scaling form proposed by Des Cloizeaux (^^): 

an = c/Nf (25a) 

10 



which under a scahng hypothesis leads to the well-known behaviour: 



pn = c3 (256) 

in dimension d = 3 

V. CORRECTIONS TO THE FLORY THEORY 

One great advantage of this method is that it allows to compute corrections to the Flory theory. 
Indeed, using eq.(9) to second order, we have: 

Z^iQ) ~ ZoiQ) exp (^-i J dr ct> (r) pQ (r) - ^ J drdr (r) Gq (r, r')(l>{r ')) (26a) 

where Gq (r, r ') is the connected Debye function^^^^ of a Brownian chain with constrained square 
radius of gyration Q : 

GQ{r,r')= / dsds'{{6{r-r{s)) 5{r' -r{s')))-{d{r-f{s))){d{r' -r{s')))) (266) 
Jo 

and is given in momentum space by (see appendix): 

Gq (k, k'^ = ^ j^J du exp (^-^(^' + ^")) (^exp (^-{^ cosu-^ smu)kk'^ - (26c) 

The parameter A in (26a) is just a remainder to keep trace of the orders of the expansion. 
Integration over the 4> variables yields: 

Z =l°°dQ exp (- (-(d- 1) In (|) + 27t'^^ - ^Tr ln((5(r- r') + Xv Gq (r,r ')) 
-| J dfdf'pQ {i^{6{r-r') + XvGQ {f,f')y'pQ (f')) 



(27a) 



Expanding to lowest order in A, we obtain: 

,2 



PF, =(3FF + x(^-j dr Gq (f, j dfdf ' pQ (f) Gq (f, f ') pq (f ')) (276) 



The first correction term is thus of the form: 



= £d„exp(-fP)(exp((2e„s„-i^sn,#) -l) (28a) 



Ar2 
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where A and B are finite numerical constants depending on the monomer size a and space dimension 
d, and the natural cut-off ^ has been introduced to avoid ultra-violet divergences. The term linear 
in N is the expected extensive part of the free energy, which in the Flory expansion appears thus as 
a correction term, and the second term is a correction to the standard Flory term in the repulsive 
energy. 

The second correction term is given by: 

J dfdf 'pQ (r) Gq (r, r ') pg [r ') = C (^^) (286) 

where C is a finite constant, depending only on the monomer size a and the dimension d . This 
term scales like the square of the Flory repulsive energy. 

More generally, it can be proven that the expansion of the free energy around the Flory theory 
will generate extensive terms, as well as powers of . This is in contrast with the usual Fixman 
expansion (^^^ , which also contains extensive terms, but which is done in powers of iV^~'^/^. The 
Flory expansion is in terms of N = A/'d+2 . 

As is often the case when calculating corrections to mean field theories, correction terms are 
much larger than the mean field contribution, in the critical region [N — > +oo). 

This allows the definition of a Ginzburg region^^*^ i.e. a typical size A^^^ax such that for 
N < iVmax, the correction terms are small compared to the Flory free energy. The precise value of 
-/Vmax depends on all the parameters of the problem, and we have not computed it explicitly, since 
it is in fact a crossover size, and its precise value is not very illuminating. 

However, it is clear that the Flory expansion will have a much larger domain of validity than 

4— d 4— d 

the Fixman expansion, since for any dimension smaller than 4, we have N'^+^ « , and 

therefore, this slower divergence of the Flory expansion might be the clue to its success. 

VI. CONCLUSION 

We have shown how the original Flory theory of polymers can be derived rigorously from 
a cumulant expansion. We obtain a Flory like free energy, with the correct original logarithmic 
term. This method can be generalized to other types of monomer interactions (e.g. bad solvent, 
polyelectrolytes, etc..) and can be applied to solutions of polymers in a straightforward manner. 
It can also certainly be useful to other classes of problems (membranes, interfaces, etc.). Finally, 
we show how this method can generate a systematic expansion around the Flory theory. The 
calculations are somewhat cumbersome, but they show that the leading corrections diverge when 
N — > +00, but much less rapidly than the usual Fixman expansion. 
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APPENDIX 

In this section, we show how to compute partition and correlation functions for a Brownian 
chain with constrained center of mass, and constrained radius of gyration. Consider the partition 
function: 

= I„ ' l" '^'') " - ^ l" '''^') (-11" ^"''*) 

(^1) 

In order to compute it, we expand the trajectories r(s) as Fourier series: 

+00 

r(s) = ^ e*T^V„ 

(^2) 



n— — 00 

1 r 

NJo 



dse ''^'^^f{s) 



The center of mass constraint impUcs that the Fourier component fo vanishes. The 6— function 
constraining the radius of gyration can be represented by its Fourier integral, and then the remaining 
Gaussian integral on r„ can be performed. After some simple algebra, we obtain: 

Zo(«=r°°|e"«n(— i^)' (A3) 
This integral can be computed by the method of residues. The poles of the function are given by: 



TV 

and they are of order d . 

Using the analytic expression for the infinite product ^^^^ : 



n°° x"^ . sinhx , 
(1 + ;^) = — W 

n=l 



we can write: 

f + CX3 



-00 

with 



f{z) = e"" ' ^ (A7) 



sinh 1 , 2 



where we have defined Q = In the following, we will forget the N factor appearing in the 
denominator of (A6), since it is just a normalization constant. We have thus: 

+ CXD 

Zo(Q) =i^res(/,z„) (^8) 

n=l 
14 



where res denotes the residue of the function. 

The exponential factor in the function / imphes terms of the form exp(— 2n^7r^g) in the sum 
(A8). 

In the case of a swollen chain, the exponent v is larger than 1/2, and thus, q — > +00 when 
N — > +00. Therefore, in the sum (A8), only the pole n = 1 will contribute, all other poles being 
exponentially subdominant. The sum (A8) reduces to: 

In the case of a collapsed chain, the exponent v is smaller than 1/2, and thus, q — > when 
N — > +00. Therefore, in the sum (AS), all poles contribute equally, but the calculation of the 
residue is simplified by the fact that q is small. Since the calculation of the residue of a pole of 
order d involves the calculation of a derivative of order d, we will specialize to dimension d = 3. A 
simple calculation of the residues yield the asymptotic formula: 

+00 

Zo(Q) = - ^(-I)"6n27r2 e'^"''^'? (^10) 

n=l 

for — > 0. Equation (AlO) can be recast in a form which emphasizes its ressemblance to the 
Jacobi elliptic Theta functions^^^^ : 

Zq{Q) = 3— ^ exp(-2nV2g + rnvr) {All) 

^ n=l 

Using the Poisson summation formula, we obtain: 

3 J +00 „ 

n — — 



n=— 00 



In the limit jj — > 0, this reduces to: 

Z,{Q) = 4 (A13) 

We now turn to the calculation of correlation functions. 

Consider the generating function for correlation functions, defined by: 

(I \ i ^ 

- - y ds P{s)j + — J ds k(s) f{s) (^14) 
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The various usual correlation functions can be obtained from (A14), by taking the proper sum of 
(5-functions for the function k{s) ( see below, eq.(A19), (A20) ) . Defining the Fourier component 
of k{s) by: 



'0 

the generating function can be rewritten as: 



kn= i ds e*T^* k{s) (A15) 
Jo 



where /(z) is given by (A7). This integral can again be evaluated by the method of residues. The 
poles arc the same as before, given by (A4), but now, due to the new exponential term, they are 
of infinite order. 

In the swollen case, q — > +C!0 , again only the pole n = 1 contributes. The calculation of the 
residue can be done and yields: 

which can be written in terms of Bessel functions*^^^-* as: 



G{kl, K) = {d- 1)1 (^^-^ j Jd-i [^2Q\k^ ^ (^18) 

The monomer density p{k) is obtained by taking: 

k{s) = N S{s- So) k {A19) 
and the pair correlation functions used in eq.(19) and (26) are obtained by taking: 

k{s) = N {6{s -so)k + 6{s - Sq) k') (^20) 
In the large distance limit, Q \ ki ^ — > 0, we obtain : 

pik) = N e"^^' {A21) 

and 

G(fc;,fcn) = e ^' 3^^"=2;r53T (^22) 

which is the form used in (19) and (26). 

In the collapsed case, q — > , all poles contribute, but in the large distance limit, the 
calculations simplify, and yield, for the density, exactly the same result as above. 
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